arXiv:1508.05349vl [physics.optics] 21 Aug 2015 


Threshold Characteristics of Slow-Light Photonic Crystal Lasers 

Weiqi XueQ Yi Yu, Luisa Ottaviano, Yaohui Chen, Elizaveta Semenova, Kresten Yvind, and Jesper Mork 

DTU Fotonik, Department of Photonics Engineering, 

Technical University of Denmark, Building SfS, DK-2800 Kongens Lyngby, Denmark 

(Dated: August 24, 2015) 

The threshold properties of photonic crystal quantum dot lasers operating in the slow-light regime 
are investigated experimentally and theoretically. Measurements show that, in contrast to conven¬ 
tional lasers, the threshold gain attains a minimum value for a specific cavity length. The experi¬ 
mental results are explained by an analytical theory for the laser threshold that takes into account 
the effects of slow light and random disorder due to unavoidable fabrication imperfections. Longer 
lasers are found to operate deeper into the slow-light region, leading to a trade-off between slow-light 
induced reduction of the mirror loss and slow-light enhancement of disorder-induced losses. 
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Slow light in photonic crystal (PhC) line-defect waveg¬ 
uides |l] enhances the interaction between the propagat¬ 
ing light wave and the material of the waveguide, and has 
enabled the demonstration of increased material nonlin¬ 
earity [2], enhanced spontaneous emission into the propa¬ 
gating mode [3, 4], and enhanced material gain .5]. Such 
engineering of fundamental materials properties is impor¬ 
tant for the development of integrated photonic circuits, 
with applications in classical as well as quantum infor¬ 
mation technology. Microcavity lasers can be realized in 
the same PhC membrane structure by exploiting high- 
quality point-defect cavities and in the past decade sig¬ 
nificant progress was made culminating in recent 

demonstrations of high-speed electrically pumped struc¬ 
tures i. Such PhC lasers allow the exploration of new 
operation regimes, such as single emitter lasing 0 and 
ultra-high speed modulation [11]. However, while it was 
shown that slow light in combination with random spa¬ 
tial disorder leads to very rich physics 12 - 18 j; the role 
of slow light on lasers realized using defect cavities has 
apparently not been systematically investigated. For the 
case of passive point-defect cavities, it is well known that 
disorder is an important factor limiting the quality fac¬ 
tor 19U221 but the role of slow light in extended active 
cavities is not well understood. 


In this paper, we report experimental results on PhC 
quantum dot lasers with variable cavity length and show 
that these attain a minimum threshold gain for a certain 
cavity length, in stark contrast to conventional lasers, 
where the threshold gain decreases monotonically with 
cavity length. We derive a rate equation including the 
effect of slow-light propagation and show that the exper¬ 
imental observations may be explained when taking into 
account disorder-induced losses. These results show that 
disorder may lead to fundamental limitations on the per¬ 
formance of nanostructured lasers, but the results also 
demonstrate a promising platform for investigating dis¬ 
order effects in active structures, such as the competition 
between deterministic cavity modes and random modes 
formed by Anderson localization [18]. 
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FIG. 1. Qualitative illustration of PhC banddiagram with 
TE-like even waveguide mode (solid line) and positions of the 
fundamental (Ml) cavity mode for L5, L10, and L20 defect 
cavities (markers). Shaded areas indicate the light cone and 
non-guided modes of the structure. The inset shows a SEM 
image of an L5 PhC cavity laser and the simulated spatial 
mode profile of the electric held amplitude. 


The PhC cavity lasers investigated here are realized 
in a 250 nm thick InP air-embedded PhC membrane 
[23] with three layers of quantum dots with densities of 
~ 5.4 x 10 10 cm -2 [ 2 ^. The PhC structure has a lat¬ 
tice constant of a = 438 nm and an air-hole radius of 
0.25a. A so-called LTV cavity jUJ is formed by omitting 
N air-holes in a W1 defect waveguide. Structures with 
different cavity lengths, i.e., L1-L20 are fabricated and 
characterized. In Fig. [lj the red solid curve illustrates 
the dispersion of the even TE-like Wl-waveguide mode 
and the corresponding modes for L5, L10, and L20 cav¬ 
ities are represented by markers. The inset shows the 
SEM image of a fabricated L5 structure and the corre¬ 
sponding simulated mode profile of the norm of the field 
amplitude in the middle of the membrane. As the cav¬ 
ity length increases, the fundamental mode moves closer 
to the Brillouin zone (BZ) edge, where the group index 
strongly increases [l], thus allowing a systematic investi¬ 
gation of the role of slow light on the laser properties. 

The samples are vertically pumped using a micro¬ 
photoluminescence set-up, with precise control of pump 
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FIG. 2. (a) Measured power of the fundamental mode Ml 
versus pump power and corresponding theoretical fit (solid 
line) for an L9 laser. The inset shows the output spectrum 
far above threshold, (b) Simulated L9 far-field profile with 
circles indicating the collection area for different NAs, and 
(c) measured emission profile corresponding to black circle. 

position and area, monitored by an infra-red camera. 
The laser emission is collected vertically by the same ob¬ 
jective, with an NA of 0.65, and monitored by a spectrum 
analyzer. The transmission efficiency from the objective 
to the spectrum analyser is ^ 30%. All measurements 
are performed at room temperature. 

Fig. Efa) shows a typical example of a measured light- 
out vs light-in characteristic, obtained by pumping an L9 
cavity with the CW output from a 1480 nm diode. The 
pump power is estimated by multiplying the excitation 
power ( P eX c ) with the ratio (Y exc ) of the PhC cavity area 
to the total pump area. The limited absorption efficiency 
of the QD layers was not taken into account. The inset 
of Fig. [2ja) displays the emission spectrum, showing the 
dominance of the fundamental lasing mode, denoted Ml, 
but also the presence of two higher order longitudinal 
cavity modes, M2 and M3. The side mode suppression 
increaes by further increase of the pump power. The solid 
line in Fig. [2^a) shows a fit to a conventional laser rate 
equation model, see e.g. [91], leading to a laser thresh¬ 
old of ^ 90 /iW and a spontaneous emission factor of 
/3 = 0.05. This /3-factor is about an order of magnitude 
smaller than values reported in [9], but enables precise 
extraction of the threshold pump power. The calculated 
far-field profile shown in Fig. [2j b) qualitatively agrees 
with measurements reported in [26J as well as our own 
experimental observations, Fig. EJc) . The calculation in¬ 
dicates a vertical collection efficiency of ^ 20% for an 
L9 cavity, which may vary slightly for different cavity 
structures. 

Fig. [3] summarizes the measured variation of the laser 
threshold with cavity length for three sample bars with 
the same PhC design. The absolute number quoted for 
the threshold pump density is defined as P exc Y exc /L , 
where L is the cavity length. While bars 1 and 2 were 
located close to each other on the wafer and fabricated in 
the same run, bar 3 is more distant and is from another 
processing batch. Gain material variations and fabrica¬ 
tion uncertainties are believed to account for the higher 


FIG. 3. Measured threshold pump power density versus cavity 
length. Three sample bars ( red “A”, blue and black “o”) 
with the same PhC design but from different positions on the 
wafer are investigated. The solid lines are guides to the eye. 


threshold of sample bar 3. The error bar indicates the 
largest variation of measured threshold between two re¬ 
peated measurements. Lasing was not observed for LI, 
L2 and L3 cavities. Despite the difference in absolute 
gain level, all the sample bars show the same qualita¬ 
tive behaviour. From L4 to L8/L9, the threshold de¬ 
creases monotonically with cavity length after which it 
increases. An optimum cavity length thus exists, cor¬ 
responding to a cavity with 8 or 9 missing holes, for 
which the laser threshold density attains a minimum. 
This behavior is strikingly different from the case of con¬ 
ventional edge emitting lasers, where the threshold den¬ 
sity decreases monotonously with length since the mirror 
losses decrease in inverse proportion to the cavity length. 

In the measurements reported in Fig. [3j we have taken 
a number of precautions to eliminate other causes for 
threshold variation. In particular, the pump profile is 
kept fixed with a diameter of 20 fim 46a), which is 
much larger than the extent of all the investigated cavi¬ 
ties from LI to L15. This ensures a uniform illumination, 
avoiding a systematic reduction of the pump efficiency 
with cavity length. Heating induced by the optical pump 
beam leads to a red-shift of the laser wavelength. The 
use of an InP membrane already significantly improves 
the thermal properties compared to more conventional 
InGaAsP structures, due to the ten times larger thermal 
conductivity of InP, [23]. To further reduce thermal ef¬ 
fect, we employed pulsed pumping with a pulse width of 
500 ns and a low duty cycle of 2%. 

The variation of the experimentally observed cavity 
mode frequencies with cavity length is shown by the 
“square” markers in Fig. |U[a). As the cavity length 
increases, all modes shift to lower frequencies and the 
spacing between neighbouring modes decreases. Quali¬ 
tatively similar behaviour was observed in Q- From the 
mode spacing and the known cavity length we can calcu¬ 
late the effective group refractive index. The extracted 
group index for the fundamental Ml mode, Fig.|U[b), in¬ 
creases monotonously with cavity length, consistent with 
the mode approaching the BZ-edge. For the optimum 
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FIG. 4. (a) Measured cavity mode frequencies (“square” 

markers) versus cavity length. The solid lines are predictions 
by the analytical theory Eq. ©• The “cross” markers are 
from the full numerical calculations. The fundamental mode, 
denoted Ml (red), is the dominating mode in all cases, and its 
threshold varies as shown in Fig 01 (b) The measured group 
refractive index (markers) for the fundamental mode Ml. The 
solid line is the calculation from Eq. ©• 


cavity length, the group index is seen to be of the order 
of - 20. 

From experimental and theoretical studies of light 
propagation in PhC waveguides, it is already known that 
disorder leads to propagation losses that strongly increase 
in the slow-light region [il FI. Since the fundamental 
mode moves deeper into the slow-light region upon in¬ 
creasing the cavity length, it is perhaps not surprising 
that an optimum cavity length exists, as demonstrated 
in Fig. 01 However, since the effective material gain per 
unit length also increases in the slow-light region [5], and 
the positive role of slow light on promotin g la ser oscil¬ 
lation in similar-type lasers was suggested 27], a more 
complete analysis is required. 

We model the laser cavity as an effective Fabry-Perot 
resonator, but with the propagating wave to be a Bloch 
mode of the W1 waveguide. This was shown in 0] to be 
a good approximation even for defect cavities as short as 
a few lattice periods. The effective length of the cavity 
is L = (TV + l)a, where N is the number of missing 
holes. By requiring that the Bloch mode-experiences a 
roundtrip phase change of an integer times 2tt during a 
cavity roundtrip [26], we find that the M'th cavity mode 
is displaced by A k = Mn/L from the BZ-edge located 
at ko = 7r /a [25[. Close to this band-minimum we may 
expand the dispersion to second-order: 


/7 . / 7T \ 1 d 2 UJ , A 7 x9 

-“(;) + 2ne [Ak) ■ (1) 

The deviation of the mode frequency from the-BZ edge, 


A/, and the corresponding light slow-down factor, S = 
n p g hc /n h g = v g /v phc , thus become 


A / 


v b gh M 2 L 
4 L 2 ’ Ml 2 ' 


( 2 ) 


Here v g (n g ) is the group velocity (group refractive in¬ 
dex) in the PhC cavity (x = phc ) and in the correspond¬ 
ing homogeneous material (x = b). We also introduced 
the characteristic length scale I 2 = ( 7r/v b g )(d 2 uj/dk 2 ), 
which is inversely proportional to an effective photon 
mass. 

Solid lines in Fig. 0] show the frequencies for the first 
five cavity modes and the group index for the fundamen¬ 
tal mode obtained from Eq. © using n b g = 3.17, I 2 = 
510 nm and a band-edge frequency of 190 THz. The 
good agreement with experiment supports the validity 
of our approximate analytical model. Mode frequencies 
(“cross” markers in Fig. 0(a)) obtained from 3D finite- 
element simulations without adjustable dispersion pa¬ 
rameters further confirm the model. 

A rate equation for the photon density including the 
effect of slow light is derived by considering the incre¬ 
mental change of the cavity photon density, N pi during 
one cavity roundtrip time, tr = 2 SL/v b . The modal 
gain per unit length of a Bloch wave propagating in the 
structure is = YSgmat , where # mat is the material 
gain coefficient for the homogeneous active medium and 
T is the confinement factor [5|. Neglecting spontaneous 
emission we find 


= v g[^9mat ( &WG + CX-MIr) / S]N p . (3) 

Here, (Xmir = Yi[l / (R1R2)] / ( 2 L) and awG denote 
the mirror and waveguide loss. We see that the spa¬ 
tial gain enhancement is cancelled out due to the longer 
roundtrip time, which is basically a statement that the 
rate of gain per unit time is unaffected by slow-light 
[I?}. Fabrication induced disorder approximately scales 
as olwg — Sql\-\- S 2 ct 2 HM3 , where the first term, with 
coefficient oq, accounts for out-of-plane scattering and 
absorption losses and the second term, with coefficient 
0 : 2 , represents back-scattering. Since the latter compo¬ 
nent does not add coherently to the laser field we assume 
that it too acts as a loss term for the photon density of 
the considered laser mode. From the rate equation for 
the laser intensity we then find that the threshold gain 
becomes 


r» = a 1+t « 2 S+^ln( :s ij-) (4) 

We see that the mirror loss is reduced in inverse 
proportion to the slow-down factor 5, while disorder- 
induced losses originating from back-scatter increase lin¬ 
early with the slow-down factor. Since S itself scales 
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FIG. 5. Calculated threshold gain, Tgth — oti versus cavity 
length for different values of 0 . 2 . I 2 — 510 nm, R\ = R 2 = 
0.98. 

linearly with L, short lasers are dominated by mirror 
loss, while long lasers are governed by disorder-induced 
losses. Fig. [5] shows the calculated length dependence of 
the modal threshold gain, Tgth — aq, for the fundamental 
mode. For 0^2 = 0 cm -1 , the threshold gain decreases 
monotonously with L, as for a conventional semiconduc¬ 
tor laser, but for non-zero 0^2, an optimum cavity length 
is found, in good qualitative agreement with the experi¬ 
mental observations. For the calculations presented here, 
we used wavelength independent mirror reflectivities of 
R \ = i?2 = 0.98, but qualitatively similar results are 
obtained using a wavelength dependent reflectivity, as 
found, e.g., in [29]. The value of will depend on the 
actual level of fabrication disorder. For our structure, 
the standard deviation of the variation of the hole radius 
was measured to be of the order of 1 nm. 

From these results we infer that, in the presence of dis¬ 
order, the quality factor of a PhC cavity is maximized for 
a specific cavity length. In order to further corroborate 
our finding we have addressed the problem using a differ¬ 
ent theoretical approach. In [30] it was shown that the 
loss of disordered waveguides can be efficiently computed 
by solving a generalized eigenvalue problem formulated in 
a basis of Bloch modes of the regular structure. Disorder 
leads to scattering among the Bloch modes and induces 
extrinsic losses of the guided modes, since these acquire 
Bloch mode components within the light cone. We ap¬ 
ply the approach of j30] to the case of a cavity, which is 
represented by a stepwise variation of the average permit¬ 
tivity along the propagation direction of the cavity, cor¬ 
responding to an effective-index approach. Considering 
only a single band, the (generalized) eigenvalue problem 
can be formulated as 

Dv/3 = UJ/3 (1 +V) V0. (5) 

Here, D is a diagonal matrix in the basis of Bloch 
modes, indexed by fc, with y/Dkk = being the waveg¬ 
uide dispersion. Furthermore, I is the identity matrix 
and V = V cavity + V dlsorder is the scattering matrix 
coupling (hybridizing) the waveguide Bloch modes due 
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FIG. 6. Calculated cavity mode spectra for L5 (left column) 
and L20 (right column) cavities in the absence of disorder 
(upper row) and with the inclusion of disorder (lower row). 
The Bloch mode distributions are shown in color-scale for 
the calculated mode eigenfrequencies (y-axis), with red (blue) 
denoting high (low) intensity. The white area indicates the 
ligt cone. 


to a determinsitic cavity “potential” as well as a ran¬ 
dom “disorder” potential. The solution of Eq. © gives 
the eigenfrequencies of the cavity modes and the corre¬ 
sponding decomposition on Bloch modes. Fig. [6] shows 
the calculated mode spectra for L5 and L20 cavities with 
and without the inclusion of disorder. A similar graph¬ 
ical representation was used in [31]. In the absence of 
disorder, the fraction of the fundamental mode that is 
within the light cone decreases with cavity length, in¬ 
dicating that the cavity loss rate decreases [25]. Upon 
including disorder, further scattering to the light cone 
takes place, reflecting the additional spatial frequency 
components introduced by the spatial disorder. 

Fig. 0 shows the relative fraction of the mode spec¬ 
trum within the light cone as a function of cavity length. 
An ensemble average over 500 realizations of the random 
scattering matrix v dlsorder was performed. We take this 
fraction as an approximate measure of the loss of that 
cavity mode. Quantitative results for the absolute loss 
rates can be obtained using the approach of [30|, but 
is beyond the scope of this paper. While the results of 
course depend on the strength of the disorder potential 
and the effective index representation of the defect-cavity, 
the features displayed by Fig. [7] are general: In the pres¬ 
ence of disorder, the light cone intensity increases, lend¬ 
ing support to the assumption of our analytical model 
that disorder-induced back-scattering increases the cav¬ 
ity mode loss. Furthermore, the mode Ml is found al¬ 
ways to be dominating. This is in accordance with our 
experimental findings, while the analytical result Eq. (|4]) 
suggests that higher-order modes become dominant deep 
into the slow-light region. 
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FIG. 7. Calculated intensity within light cone versus cav¬ 
ity length for the fundamental mode Ml, with and without 
disorder, and for the higher-order mode M2, with disorder. 


In conclusion, we have experimentally shown that pho¬ 
tonic crystal lasers display intriguing threshold character¬ 
istics, with a minimum threshold gain being attained at 
a specific cavity length. An analytical model was derived 
that well accounts for the experimental results, suggest¬ 
ing that the observed threshold behavior is a result of the 
interplay between slow light and disorder-induced losses. 
These results are important for the understanding and 
further development of photonic crystal lasers. The re¬ 
sults also demonstrate a promising platform for system¬ 
atic experimental and theoretical investigations of disor¬ 
der effects in active (amplifying) media and in cavities, 
which remain relatively unexplored compared to passive 
photonic crystal waveguides. For instance, a controlled 
amount of intentional disorder may be added lji, 21] and 


the relative strength of deterministic cavity modes ver¬ 
sus random Anderson-localized modes [lj| might be con¬ 
trolled via the mirror reflectivities. 
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